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The behavior of the renormalized spin 2-point function in the 0(3) and dodecahedron spin model 
are investigated numerically. The Monte Carlo data show excellent agreement between the two 
models. The short distance behavior comes very close to standard theoretical expectations, yet it 
differs significantly from it. A possible explanation of this situation is offered. 
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Are two dimensional (2D) nonabelian cr-models asymp- 
totically free? After more than two decades, the answer 
remains unclear, with lots of supportive evidence on both 
sides of the issue. The main point of this letter is to pro- 
vide some fresh evidence that even though the orthodox 
picture provides an excellent description of physics at in- 
termediate momemta and distances not too short, it is 
in fact incorrect, a situation which I belive occurs also in 
QCDi. 

To that end I will present the result of Monte Carlo 
(MC) investigations of the continuum limit of the spin 
2-point function in four spin models: i) Ising model, ii) 
0(2) model, iii) 0(3) model and iv) dodecahedron spin 
model. As I will show, there is strong numerical evi- 
dence that the continuum limit of the dodecahedron spin 
model is identical to that of the 0(3) model. Under- 
going a freezing transition at nonzero temperature, the 
latter is not likely to be asymptotically free. Moreover 
the continuum limit of the spin 2-point function of the 
0(3) model disagrees with the PT formula with the value 
of A/m predicted by Hasenfratz, Maggiore and Nieder- 
mayer based on the thermodynamic Bethe ansatz [0 and 
the normalization constant fixed by the scaling hypoth- 
esis of the form factor approach proposed by Balog and 
Niedermaier (BN) [|. 

The numerical data suggest that in all the four models 
investigated, the short distance behaviour of the spin 2- 
point function is of the form 

As a biproduct of my investigation, I report also on 
the behavior of the critical exponent ry, which probably 
equals 1/4 in all four models. The subject of logarith- 
mic corrections to = 1/4 in the 0(2) model has been 
discussed in several recent papers (3|, [D, since 
a theoretical prediction stemming from the perturbative 
renormalization group approach exists and it is not 
clear that the data are consistent with it. As I will dis- 
cuss shortly, the existence of a continuum limit obtain- 
able from the lattice model via multiplicative renormal- 
ization, an assumption which seems to be corroborated 
by the numerics, suggests that no such logarithmic cor- 
rections are present. 

I would like to begin by reviewing the intimate con- 



nection between the short distance behavior of the con- 
tinuum spin 2-point function and the critical exponent rf. 
First some definitions: let P — {p, 0) and 
- Spin 2-point function G{p): 



G{p) = ^{\s{p)n §{p) 



Magnetic susceptibility x- 

X = G{0) 

Correlation length ^: 
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As the continuum spin 2-point function I will use the 
limit ^ — > oo of 



G{x/0 ^ {s{0).s{x))e/x 



(4) 



Let z be the physical continuum distance {z = x/^, limit 
^ oo); the exponents fj and r are defined by requiring 
that in the limit z ^ the expression 



Giz)z'>{logiz)Y 



(5) 



is finite and different from 0. Let rj and r be defined so 
that in the limit ^ — > oo the quantity 



Ckt 



X 



(6) 



is finite and different from 0. 

I would like to claim that the existence of a continuum 
limit with the short distance behaviour described by (5) 
requires rj — fj and r — —f. Indeed in the continuum 
limit, the limit z — > is reached by letting both x oo 
and ^ — > oo but x/(, —^ 0. In fact the same considera- 
tion shows that in the same limit the lattice spin 2-point 
function must have the property 



{s{0).s{x))x'' 



(7) 
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goes to a finite nonzero constant (i.e. there are no loga- 
rithmic modifications). 

The Kosterhtz prediction is that f — —1/16. Fi- 
nite size scaling at the fctpoint do indicate a negative r 
though smaller in magnitude . The paradox discussed 
by many § , 1^ , |6| is that the thermodynamic data sug- 
gest a positive value for r |4|, p. As explained above 
this is consistent with the fact that r — —f, rather than 
r = f as incorrectly claimed in the literature. 

However the fact that eq.(7) requires the lattice spin 
2-point function to decay as a pure power with no loga- 
rithmic corrections suggests that most likely r = — f = 0. 
Indeed while for ^ < oo eq.(7) is valid only for a; ^ ^, 
baring some nonuniformity, eq.(7) ought to be describing 
the large distance behaviour of the lattice spin 2-point 
function at the critical point. 

Next I would like to discuss the main point of this pa- 
per, which concerns the asymptotic behaviour of the con- 
tinuum spin 2-point function G{z) at short distances. It 
is known that in the Ising model i] — 1/4 r — 0. I stated 
above what perturbative renormalization group (RG) ar- 
guments predict for 0(2). For 0(3) similar arguments ||] 
predict ry = r = — 1. However when combined with the 
thermodynamic Bethe ansatz and the scaling hypothesis 
in the form factor approach, this latter prediction be- 
comes even more precise : 



G{z) 



1.000034657 

37r3 1.0016872 
1.1159 + ln(t)/t 



it 



ln{t) 

.1159/i)^ 
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where t = —ln{zK),K = e/8. I am not aware of any 
predictions regarding rj and r for the dodecahedron spin 
model. 

In Fig.|l| I present my Monte Carlo results for G{z)z^/'^ 
for (in rising order) Ising, 0(3) and 0(2); I also show the 
curve for free field (lowest curve) and the PT curve with 
the normalization constant predicted by Balog and Nie- 
dermaier (upper curve) The MC data were taken at 
^ ^ 167 and L=1230. They represent the on axis corre- 
lation and I will discuss shortly the lattice artefacts. For 
now I would like to make the following observations: 

1. In all three models G{z) approaches the free field 
value for z —f oo, as required by the Orstein-Zernke be- 
haviour 1^. 

2. For z > .05 G{z) in the three models behaves quite 
similarly. 

3. The data suggest that in all three models rj = 1/4, 
but that compared to the Ising model the short distance 
behavior may be more singular in 0(2) and less singular 
in 0(3). This is would be consistent with a negative f in 
0(2) but a positive one in 0(3) (although the data are 
also consistent with f = 0) . 

4. Even though for .02 < z < .1 the difference with the 
BN refined PT curve is small (about 2%), it is clearly 
there and, as I will indicate next, it is not a lattice arte- 
fact. 
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FIG. 1. Spin 2-point function G{z) versus physical distance 
z for the Ising, 0(3) and 0{2) models. Also shown the free 
field curve and the BN refined PT curve for 0(3). 
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FIG. 2. Spin 2-point function G{z) versus physical distance 
z for the 0(3) model at several correlation lengths. Also 
shown the BN refined PT curve. 

In Fig.^ I present the MC results for the 0(3) model 
for ^ ranging from 65 to 310. The error bars are about 
.3%, too small to show. The data were produced using 
Wolff's method jl^]. At each /3 and L value, at least 
five independent runs were made and the error estimated 
using the jack-knife method. The primary source of error 
are the values of x ^-nd ^, whose determination involves 
large distances. To give a feeling for the size of the lattice 
artefacts, for ^ ^ 310 the (solid) curve represents G{z) 
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on axis while the points G{z) along the diagonal. Both 
this test and the good agreement of these data with the 
data at ^ ~ 230 suggest that, barring some very slow 
drift, the continuum limit has been reached for z > .03. 
The continuum values are clearly off the BN refined PT 
curve (upper curve). 

As explained above, the critical behaviour of x ^^^id ^ 
carry relevant information regarding the short distance 
behaviour of G(z). In Fig.^I represent ln{x) — l-75ln{S,) . 
The upper curve represents the MC data. The lower 
(dashed) curve represents the ratio of the values of this 
quantity computed via the BN formula 



37r 1.0016872 
41.000034657/32 



.1816 .133 



.1362 



/32 



/33 } 
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over its MC values. A correct prediction of the BN for- 
mula would require the dashed curve to approach 1 for 
^ ^ cxo. The data produce a (dashed) curve which is not 
flat and which seems to approache the line at 1 with a 
nonzero slope. 
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FIG. 3. ln{x) - 1.75Zn(0 versus ln{^) for the 0(3) model. 
Dashed curve represents MC data combined with the BN 
bootstrap prediction, a correct prediction requiring conver- 
gence to 1 for ^ oo. 

The upper (solid) curve indicates that r could be dif- 
ferent from and in fact negative. To determine whether 
there is a logarithmic modification to 77 = 1/4 and ob- 
tain a precise value for r one would need data at larger 
^. However the lack of smoothness of the curve suggests 
that the errors are underestimated already at ^ ^ 310 
(probably due to critical slowing down). 

The last piece of numerical evidence which suggests 
that 0(3) is not asymptotically free comes from a com- 
parison of its spin 2-point function with that of of the 



dodecahedron spin model. This is shown in Fig. which 
shows G{z) for the latter model at ^ ~ 11, 19, 34, 66 and 
120 (broken hne) and G{z) for the 0(3) model at ^ - 122 
(solid line on axis, dots on the diagonal). As the correla- 
tion length increases, G{z) for the dodecahedron model 
approaches that of 0(3) and at ^ ^ 120 the two curves 
are practically indistinguishable for z > .04. It appears 
rather certain that the two models share the same G{z) 
in the continuum limit. However, as already mentioned, 
it is highly unlikely that the dodecahedron model, under- 
going freezing at finite (3, exhibits asymptotic freedom. 
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FIG. 4. Spin 2-point function G{z) in the dodecahedron 
model for several correlation lengths. Also shown G{z) for 
the 0(3) model at ^ - 122 and the BN refined PT curve. 

Discussion 

1. I believe the numerical evidence presented above 
suggests that in the 0(3) nonlinear a model the short 
distance behaviour of G{z) corresponds to 77 = 1/4, f > 
0. 

2. Perturbative renormalization group arguments pre- 
dict that in all 0{N) models N > 3 r] — 0, f is negative 
and varies monotonically with N from -1 for 0(3) to - 
1/2 for O(cxd). Based on the observed behavior of 0(3) I 
would conjecture that for all 3 < /V < 00 77 = 1/4. Since 
the spherical model O(oo) has 77 = 0, r ^1/2, such 
a scenario would correspond to a nonuniformity of the 
1/iV expansion. The uniformity of the 1/N expansion 
for /? — > 00 has been questioned before |ll| . 

3. Since for .02 < z < .1 the deviation from the BN 
refined PT prediction is 2% or less, one may wonder why 
that is so? I think the answer is in Fig.l, where I com- 
pare G{z) for different models. Please note that the free 
field curve describes also the spherical model and G{z) 
for it is quite close to that of 0(3) for z not too small. 
It is well known |l^ that the 1/N expansion is legiti- 
mate at fixed P = jj: and the only issue jl^ is whether 
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this expansion is uniform for P —>■ oo. The reason the 
orthodoxy comes so cfose to the truth is two fold: i) the 
fast convergence of the 1/N expansion for moderate /3 
and ii) the large value of /3crt- Thus, even though I be- 
lieve both the expansion and perturbation theory 
are nonuniform, they come very close to the truth at 
moderate correlation length. Since for distances not too 
small or momenta not too large the continuum limit is 
well described at moderate ^ values, the orthodoxy will 
describe correctly physics at such distances and/or mo- 
menta. It is only at short distances or large momenta 
that the true test of the orthodoxy can be performed. 

4. Until now I discussed only some new numerical 
information. Some readers, such as those who recently 
wrote 'We believe however that the continuing accumu- 
lation of unambiguous, consistent and increasingly accu- 
rate numerical support for the RG predictions from a va- 
riety of independent approaches leaves liitle if any space 
for alternative pictures. ' p^ , may not find my numerical 
evidence conclusive. I will invite them then to reconsider 
the two theoretical arguments advanced by Patrascioiu 
and Seller regarding the existence of a massless phase in 
all 0{N) models: 

i) We showed rigorously that if a certain cluster, 
baptized the equatorial cluster, does not percolate, the 
model must be massless. We also advanced some nonrig- 
orous arguments why that cluster did not seem likely to 
percolate After seven years since those arguments 
were advanced, no mathematical physicist has proved the 
contrary or given us an example of how this equatorial 
cluster might percolate (see Open Problems in Mathemat- 
ical Physics, http://www.iamp.org). Can any skeptical 
reader give such an example? 

ii) We also showed rigorously that among smooth con- 
figurations, the dominant ones must be the superinstan- 
tons []l6| not the much publicized instantons. But in a 
superinstanton configuration there is a ring of arbitrary 
size in which the spin points roughly in a certain direc- 
tion. Hence in such a configuration, the inverse image of 
a small patch of the sphere will form a ring and hence the 
equatorial cluster will not percolate. Therefore the exis- 
tence of superinstantons as the dominant configuration at 
weak coupling reinforces the argument for the existence 
of a massless phase and hence the absence of asymptotic 
freedom in the massive phase. 

I think irrespective of any numerics, to paraphrase 
Butera and Comi |^3|, 'the two arguments mentioned 
above leave little if any space for the standard picture 
to be correct'. 

5. Even though strictly incorrect, the standard pic- 
ture provides an excellent phenomenological description 
at distances not too small {z > .02) or momenta not 
too large {p < 27r/.02 ^ 300. This is not something 
never seen before. Indeed it is taught in every course in 
Quantum Mechanics that the probability to find a parti- 
cle trapped in a potential well (say an a-particle) decays 
exponentially in time. The time constant of that expo- 



nential is called the life time of the state, and innumer- 
able experiments over the years have measured such life 
times. I am not aware of anybody ever having claimed 
to have detected deviations from this alleged exponential 
law. Yet it is a mathematical fact that if the particle 
decays at all, the large time asymptotic behavior cannot 
be exponential, but some inverse power of the time [ pT[ . 
This fact is little known in the physics community (wit- 
ness the many papers by particle physicists computing 
the lifetime of the 'false cosmological vacuum') because 
in fact the exponential law describes very well the data. 
My claim is that a similar situation occurs in the 2D 
nonlinear a models (and probably QCD4): even though 
eventually, at /3 sufficiently large superinstantons win and 
render the model massless, at moderate values of (3 such 
as /3 = 2 for 0(3), instantons (localized defects) domi- 
nate the typical configuration, which behaves very much 
as predicted by perturbation theory and/or the 1/A^ ex- 
pansion. Only by probing sufficiently small distances or 
large momenta should one be able to detect major devi- 
ations from the expected behavior. 

These ideas stem from my long term collaboration with 
Erhard Seller and were crystalized by my recent collabo- 
ration with J.Balog, M.Niedermaier, F.Niedermayer and 
P.Weisz. 
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